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APPENDIX A 

4e 4e 4t 3(e afe :|c lie * 4c 4e 4e :(c ^ 

#include <stdio.h> 
#mclude <stdlib.h> 
* ^.r^rM^p"^ #include <striiig.h> 
#mclude <stdlib.h> 
#mclude <ctype.h> 
#include < assert. h> 
#iBclude <bitvector.h> 
#inchide <list.h> 



, In 



m 



struct Procedure : public Object { 
char *code ; 
char *desc ; 
double fee ; 

int isPanel ; // indicates whether procedure is a panel 
BitVector panelProcs ; // if is panel, this contains the panel codes 

>; 

struct TreeNode : public Object { 

BitVector procedures ; // hst of procedures for this branch 
List branches ; 

}; 

int numProcedures = 0 ; 

Procedure *procedTires = NULL ; // array of procedures 
TreeNode rootNode ; 



int findCode(char *code) 
{ 

for(int i=0;i<nuniProcedures;i+-f ) 
if( strcmpi(procedures[i].code,code) ==0) return i ; 
return -1 ; 
} 

void trimString(char *str) 
{ 

char *p = str ; 
while(isspace(*p)) p++ ; 
int len = strlen(p) ; 
memcpy(str,p,len) ; 
str[len] = 0 ; 
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len = strlen(str) ; 
if ( len ) { 

p = str+(strlen(str)-l) ; 
while(p>=str &&isspace(*p)) *p~ = 0 ; 
} 
} 

void loadPanelDefinitionsCchar *filenaine) 

{ 

FILE *fp = fopen(mename,V) ; 
assert(fp) ; 

char buf[1024] ; 

// load the panel definitions 
while(fgets(bTif,sizeof(buf),fp)) { 
trimStringCbuf) ; 

// get code and procedure 
char *p = strtokCbuf,", ") ; 
int id = findCode(p) ; 
if(id>=0){ 

Procedure *panel = &procedures[id] ; 
panel- >isPanel = TRUE ; 
while(p=strtok(NXJLL;*, ")) { 

int pid = fmdCode(p) ; 

if(pid>=0) 

panel- >panelProcs.set(pid) ; 

} 
} 
} 

fclose(^) ; 
} 

void buildProcedureListCchar *filename) 
{ 

FILE *fp = fopen(filename,"r") ; 
assert(fp) ; 

char bufI1024] ; 
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// get the first line containing the number of procedures 
fgets(buf,sizeof(buf),fp) ; 
numProcedures = atoi(buf) ; 

assert(numProcedures>0&&numProcedures<20000) ; // the 20000 cap is 
arbitrary just for sanity's sake 

// allocate the array of codes 

procedures = new Procedure[numProcedures] ; 

assert(procedures) ; 

// load the procedures 
int i = 0 ; 

while(fgets(buf,sizeof[buf),fp)) { 
trimStringObuf) ; 

// get code 

char *p = strtok(buf,", ") ; 
procedures[i].code = strdup(p) ; 
procedures[i].isPanel = FALSE ; 

11 ^qX, fee 

p = strtok(NULL,", ") ; 
procedures[i].fee = atof(p) ; 

// get desc 

p = strtok(NULL,"") ; 
if(p){ 

while(*p&&(isspace(*p)| l*p==';)) p++ ; 
procedures[i].desc = strdup(p) ; 

}else 

procedures[i].desc = strdupC") ; 

//printf("%s,%f,%s\n",procedures[i] .code,procedures[i] .fee,procediires[i] .des 
0; 

i++ ; 
} 

fclose(fp) ; 
> 

BitVector *loadOrders(char *filename) 
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{ 

BitVector *orders = new BitVector ; 
assert(orders) ; 

FILE *fp = fopen(filenaine,"r") ; 
assert(fp) ; 

charbiif[1024] ; 

// load the panel definitions 
while(fgets(buf,sizeof(buf),fp)) { 
.trimStringCbuf) ; 

// get code and procedure 
char *p = strtokCbuf,", ") ; 
int id = findCode(p) ; 

assert(id> =0) ; // Invalid procedure code in ordered procedures 

if( id > = 0 ) orders- >set(id) ; 

} 

fclose(fp) ; 
return orders ; 
} 

void unbundlePanels(BitVector *procList) 
{ 

// this routine makes multiple passes through orders list in case 
// there is an occurance of a panel within a panel 
int changed ; 
do{ 

changed = FALSE ; 

for(int i=procList->start();i<=procList->end();i++) { 
if( procList->get(i) && 
procedures[i].isPanel ) { 
procList-> clear (i) ; 

procList- > or (&procedures[i] .panelProcs) ; 
changed = TRUE ; 

} 
} 

} while(changed) ; 
} 
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void dumpPrdcedureList(BitVector *procList) 
{ 

int nximPanels = 0 ; 
int nximProcs = 0 ; 
double totalFee = 0.0 ; 

// dximp ordered procedures 

forCint i=procList->start();i<=procList->end();i++) { 
if( procList->get(i) ) { 
nuiiiProcs++ ; 

totalFee += procedures[i].fee ; 

if( procedures[i]isPanel ) nuinPanels++ ; 

priiitf("%s - %s %s\n", 

procedures[i] .code,procedures[i] .desc,(char 
*)((procedures[i].isPanel)?"(PANEL)'^:'"')) ; 
} 
} 

printfCVi") ; 

printf ("Total panels: %d\n", numPanels) ; 

printfCTotal procedures: %d\n", numProcs) ; 

printfCTotal Fee for Listed Procedures: $%-5.2i\n", totalFee) ; 

} 

int panelProceduresMatch(Procedure *panelProc,BitVector *procList) 
{ 

// first check to see if panel proc is special chemistry procedure 
BitVector bv ; 

bv.or(&panelProc->panelProcs) ; 
bv.and(procList) ; 

//printf("****Dumping merged procedure list****\n\n") ; 
//printfCTanel Procs\n") ; 

//dumpProcedureListC&panelProc- >panelProcs) ; 
//printf("procList\n") ; 
//dumpProcedureList(procList) ; 
//printf("Merged\n") ; 
//dumpProcedureList(&bv) ; 
//printf("**** END OF DUMP ****\n\n") ; 
for(int i=bv.startO;i<=bv.end();i+-h) { 
if( bv.get(i) ) { 

printf(" Matched %s on procedure 
%s\n",panelProc- > code,procedures[i] .code) ; 
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return TRUE ; 

} 
} 

return FALSE ; 
} 

void dtiinpTree(TreeNode *node,int Ivl) 
{ 

TreeNode *pNode = (TreeNode *)node->branches,First() ; 
int cnt = 0 ; 
while(pNode) { 
printf("%*.*s%d:%d:",lvl,lvl," 
Mvl,cnt) ; 

for(int i=pNode->procedures.start();i<=pNode->procediires.end();i++) { 
if( pNode-> procedures. get (i) ) 
printf ("%s,",procedures[i] .code) ; 

} 

prmtf("*\n") ; 
dumpTree(pNode,lvl+l) ; 
cnt ++ ; 

pNode = (TreeNode *)pNode->next ; 

} 

} 

void buildTreeCTreeNode *node,BitVector *procList) 
{ 

TreeNode *pNode = new TreeNode ; 
assert(pNode) ; 

pNode->procedures.or(procList) ; 
node->branches.Insert(pNode) ; 
for(int i=0;i<niunProcedures;i++) { 

if[ procedures[i].isPanel && 

panelProceduresMatch(&procedures[i],procList) ) { 
// printfCMatched %s\n",procedures[i].desc) ; 

BitVector panelProcList(nuniProcedures,TRUE) ; 

panelProcList.clearAllO ; 

paneIProcList.or(&procedures[i] .panelProcs) ; 

panelProcList.and(procList) ; 

pNode = new TreeNode ; 

assert(pNode) ; 

pNode->procedures.set(i) ; // set the panel ID 
node- > branches.Insert (pNode) ; 
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Bit Vector remainingProcList(numProcedures,TRUE) ; 
remainingProcList.or(procList) ; 

forCint j=panelProcList.start();j< =panelProcList.end() J+ +) 
if( paiielProcList.get(j) ) remaiiimgProcList.clear(j) ; 

if( IremainingProcList.isNuUO ) 
buildTree(pNode,&reiiiainingProcList) ; 

} 
} 
} 

double calcFee(BitVector *procIist) 

{ 

double totalFee = 0.0 ; 

// dump ordered procedures 

forCint i=procList->start();i< =procList->end();i+ +) 

if( procList->get(i) ) 

totalFee += procedures[i].fee ; 
return totalFee ; 
} 

double getMiniinumF ee(TreeNode *node,BitVector *minProcs) 

{ 

TreeNode *pNode == (TreeNode *)node->branches.First() ; 
TreeNode *minNode = NULL ; 
double fee = 0.0 ; 

BitVector imnVector(numProcedures,TRUE) ; 
while(pNode) { 

double val = calcFee(&pNode->procedures) ; 
BitVector bv(numProcedures,TRUE) ; 
bv.clearAllO ; 

bv.or(&pNode-> procedures) ; 

val = val + getMinimiiinFee(pNode,&bv) ; 

if( IminNode 1 1 val < fee ) { 

minNode = pNode ; 

fee = val ; 

minVector.clearAllO ; 
minVector.orC&bv) ; 
} 

pNode = (TreeNode *)pNode->next ; 
} 
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minProcs->or(&ininVector) ; 

return fee ; 

} 

int main(int argc, char *argvG) 
{ 

printf ("Dynamedix Medicare Billing Optimization Patent Algorithni\nUsing 
Dallas Locale and Carrier Code For Fee Structure") ; 

buildProcedureList("codes.txt") ; 
loadPanelDefinitions("panels.txt") ; 

// load the orders, build the tree and report the results 
BitVector *orders = loadOrders("orders.txt") ; 

// Unbundle and dump 
printf('^nOriginal Order\n\n") ; 
dumpProcedureList(orders) ; 
unbxmdlePanels(orders) ; 

printf(''\nOriginal Order Unbundled with Duplicates Removed\n\n") ; 
dxmipProcedureList(orders) ; 

printf C^n\nBuilding Tree Structure\n") ; 
TreeNode root ; 
buildTree(&root,orders) ; 
dumpTree(&root,l) ; 

BitVector optOrder(numProcedures,TRUE) ; 
double fee = getMinimumFee(&root,&optOrder) ; 
printfC^n\nOptimum Order\n") ; 
dumpProcedureList(&optOrder) ; 

return 0 ; 
} 
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